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Introduction 

It is well-known, that in some cases the contributions of the high-dimensional (d > 6) 
operators in the QCD sum rules [1] are important. For estimation of the vacuum average 
of high- dimensional operators usually is used factorization hypothesis [1]. According 
to this hypothesis it is assumed that in the high-dimensional operators expansion over 
intermediate states, vacuum states contributions are dominant. Formally this assumption 
can be written as 

< Oi ■ 2 >— < Oi> ■ < 2 > 

where < O > are some color scalar operators. 

In [1] factorization hypothesis was used to estimate vacuum average of the 4-quarks 
operator (d — 6). Gluons operators with d — 6 and d — 8 was been calculated in [2,3] 
(see also discussion about it in [4]). But for operators with dimension (d — 7) a large 
number of new vacuum averages appears, and some of them can't be reduced to a product 
of vacuum averages of the operator with lower dimension. Nevertheless in this paper we 
make an attempt to estimate vacuum averages of all operators with d = 7. The method 
we offer is based on factorization hypothesis. 

It is well-known, that though factorization hypothesis is confirmed in I/TV limit [1], but 
for real world with N=3 factorization hypothesis in some cases has rather bad accuracy 
(see for example [5]). So in this work we shall use factorization hypothesis in maximal 
"soft" form, i.e. we suppose: 

1) If any operator can be saturated by vacuum intermediate states (i.e. if vacuum inter- 
mediate states contribution exist and is not zero ) we estimate it by its factorized value 
and suppose accuracy 30%. 

2) To estimate other operators one should try to express them through factorizable oper- 
ators, if it is possible. 

To avoid uncertainty, we also suppose: 

3) Vacuum average of high- dimensional operators should not depend on the way of fac- 
torization, and, particularly, vacuum average of operators, containing derivatives should 
not depend on the fact are the equation of motion taken in account before or after fac- 
torization. 

This assumption 3) can be treated as a condition of self-consistence of factorization hy- 
pothesis. This three assumptions we hereinafter shall call anzatz. 

The main idea is to consider the vacuum average of the operator with large dimension 
(d = 10). Using fagtorization hypothesis one can express it in terms of the product of 
vacuum averages of operators with lower dimension (d < 7), some of which are known and 
some are unknown. According condition of self-agreement of factorization one can write 
a number of relations for this unknown operators and estimate them. This method will 
be explicitly explained in section 1. As by-product some estimations for vacuum average 
of the operators of dimension 10 became available. Of course one must note that this is 
only some rough phenomenological estimations and do not claim for high accuracy. An 
accuracy of our estimation are about a factor 2. 

The paper is organized as follows. 
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In section 1 we describe the method on an example of calculation vacuum averages of 
the of dimension 7 operators, constructed from quark and gluons. 

In section 2 we discuss vacuum averages of the 7-dimension operators with one deriva- 
tives. So it appears possible to evaluate all vacuum averages of the dimension 7. In 
section 3 obtained results are used to calculate high dimensional operator contribution to 
Bjorken [6] and Ellis- Jaffe [7] sum rules. 
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Section 1 



In this section we'll discuss following vacuum averages of the of dimension 7 nonfac- 
torizable operators: 

< g 2 gd nkl \ l Gl v G%q > < g 2 gf nkl X l G; u G k ^q > 
Rd ~ 24 ; Rf ~ 24 

c i < g 2 qi 5 G^G^q > i < g 2 q^dr kl X l G n ^ v q > 

Si- ; d ~ 24 (1) 

where G™ v = G™ lVl ■ 

For convenience hereinafter following notations will be used 

RdJ = RdJ <qq> Ri =< g 2 G 2 >< qq > 2 /24 

Si, d = S 1>d <qq> N = (< gqG^q >) 2 /24 (la) 

_ _ r^n . /~i2 /~m /~m 

^*liv — 2 /ii" ~ ^ iiv^ \xv 

Note, that at standard choose of gluon, quark and quark-gluon condensates R ~ iV . 
We assume that in our accuracy vacuum averages for u and d-quarks are the same, so, 
for example 

< g 2 u<T kl \ l Gl v G\ v u >=< g 2 dd nkl X l G; u G k u d > 

The method of estimation R d , Rf--- is based on our anzatz. We consider some vacuum 
averages of the of dimension 10 operators and factorize them in two different ways. Once 
we use equation of motion before factorization, and other time -after. We require this 
two results to be the same (with accuracy ~ 30 %, which is the accuracy of factorization 
itself). Because of such uncertainty in order to estimate vacuum averages (1) one have 
to use only those equations, which are enough (> 50%) sensitive to unknown values of 
vacuum averages (1). 

We shall illustrate the method on an example of the following operator with dimension 

10 

< Ti >= - < (d[V 2 u}) ■ (u[V 2 d}) > 

where square of a bracket like [V 2 w] mean that a derivative acts only on the quark operator 
in bracket. From one side < 7\ > can be immediately factorized in the single way 

< Ti >= — < u[V 2 u] >< d[V 2 d] > (2) 
Using equations of motion Vg = and the fact that 

V 2 = VV + \gG^\ (where V = 7 ^) (3) 
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One can find 



< T 1 >=--N 



(4) 



where N is denoted in (la). 
From other side if take into consideration equation of mouton from very beginning, one 

get 



<T 1 >=- 9 -< d^u ■ u^d ■ G n ap G% > 



(5) 



Now one can do here Fiertz transformation (for a simplicity we will write down it only 
for color indexes, and for Lorenz indexes it is meant) 



< T x >= 9 -\ l - < (u T ,u T ) ■ (d p ,\ k \ n d p )G k a ,G; £ > + 



+ \ < (u T 'X l u T ) ■ (d p >\ k \ l \ n d p )G k /3 G™ £ > j • \ ■ af T a% 



Using the fact that 



ab 



X k \ l \ n ) = 1 ( d nkl + if nkl ) 5 ab + ^ ( (X k ) ab 5 ln + (X n ) ab 5 lk 



13 
21 



(6) 



21 



^l^ab ^kn _|_ Qmkln ^ym\ab 



(7) 



(where O mkln = Tr(X m X k X l X n )/2 - §(5 mk 5 ln + 5 mn 5 kl - ^5 ml 5 kn ) - is traceless matrixes 
by each pair of indexes) we find 



^ rp 9 a/3 <pe j 1 



< (u T ,u T ) ■ (d Pl X k X n d p G k ^ £ ) > 



+i < {u T ,X\d nkl + ir kl )G k ap Gl £ u T ) ■ (d p> d p ) > 



26 
+21 



< (u^G^Ur) ■ (d pl G a pd p ) > + < (u T iG a/ 3U T )(d p 'G ve d p ) > 
5 



42 



< {u T ,\ l u T ){d p ,\ l d p )G k Gl e > 



+ -< {u T ,\ l u T ){d p ,O mkln \ m d p )G k ap Gl £ > 



(8) 



Fourth term in (8) after factorization appears to be expressed in terms < (uX l Tu) ■ 
(dX l Td) >< g 2 G 2 > where T = 1, 

75, 7a* 5 757^5 oilv ■ This terms assumed to be 0, 
because, as was shown in [5], < [uX l Yu] ■ (dX l Td) > are negligible small for T = 1,7^ and 
there is not any reasons to believe that for other T result will grow on an order, so, we 
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can neglect by fourth term. The last, fifth, term in (8) appears to be at factorization 
due to fact that O nklm is traceless matrixes The first three terms in (8) allows vacuum 
intermediate state, so < Ti > can be factorized and after some calculations one can found 

<T 1 >=R 1 /6 + R d /2 + R f + ^N (9) 

Here we omit terms, proportional to Si, S d , because, as will be shown later, they are 
negligibly small. Comparing (4) and (9) we get 

N/2 = RJ6 + R d /2 + R f + ^N (10) 

Note that this equality is rather sensitive to unknown vacuum averages R d and R s . Really, 
if we suppose R d = Rf = 0, then left and right side differ from each other more than 
twice. 

One more equation can be obtained, if we consider such an operator 

< T 2 >=< d x [L^u] p ) ■ (u T ,u T ) > Y T ■ (i a Y)x P (11) 
where L a » = (V^V^ - V a V M V^) . 

< T 2 > can be factorized immediately as < T 2 >= - < u^f a Y[U. afl u] >< dd > ~. 
Using equations of motion after some calculations one can find 

<T 2 >=Ri/3 + R d /2 + R f /2 (12) 
From other side by use of equation of motion one can write (11) in the form: 



< T 2 >= -g 2 < {dG^G^l^iud) > 
Then after Fiertz transformation just as in previous case one can get 

<T 2 >= R^S + Rd + Rf + ^N (13) 



From (12) and (13) we have 

Ri/3 + {R f + R d )/2 = ^- + R d + R f + ^. N (14) 

Note, that (14) also is sensitive to R d , Rf. 
Using the same procedure for the operator 

< w >= - < i(d l5 e a ^[V a VpVpV v u])(ud) > (15) 
we get such an equation: 

Si/3 + S d /2 = Si/3 + S d -^N (16) 

From (16) we can conclude S d ~ iV/5 <C N so we put it 0. What about Si, we assume it 
to be negligibly small too. Really, it is easy to show by direct calculations 
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Si = 24 < qq ■ G 2 > -- < qortov G^G^q > (17) 

So Si is equal to difference of two factorizable vacuum averages which cancel each other 
after factorization. So we can expect that Si should be much more less then this vacuum 
averages (which are of order of Ri) and we can put it within our accuracy. 
From (10), (14) we found 

Rf ~ 4/21 • iV — Ri/3 <C R\ (or N); so E ; ~0 

R d ~ -26/21 ■N + R l /3~ -N (or R d R x because i? x ~ N, see (1)) 

Finally, we write down results of this section: 

R d ~ -R u ~ S d ~ « (18) 

Section 2 

In this section we will discuss vacuum averages of the four-quark operators with one 
derivative, such us: 

X, =< (g[V a g])(g 7 °g) >; *i =< (qi 5 lV a q])(q^ a q) > 
X 2 =< (gA fe [V a g])(gAV<?) >; X 2 =< (qX k i 5 [V a q])(q\ k l5l a q) > 
Yi = i < {qi T [V £ q]){qa T£ q) >; Y x =i< (g7s7 T [V e g])(g75^ T£ g) > 

Y 2 = i < (q\ k Y[V e q])(q\ k a T£ q) >; Y 2 = i < (q\ k l,l T [V £ q])(q\ k l5 a T£ q) > (19) 

Note, that one can't factorize this operators immediately (after proper Fiertz transforma- 
tion), because due to equation of motion they became zero after factorization (see point 
1,2 of our anzatz). Note also, that all other vacuum averages of the four-quark operators 
with one derivative easily can be expressed by this eight, by help of equations of motions. 
Of course, this eight operators aren't independent. First it can easily be shown, that 
Y\ = —Xi. Really, due to C-parity 

Yi = \ < {qi T [V £ q} + \qV £ ]l T q) • (qa Te q) >=\< [d £ (qYq)] • (q<r re q) > (20) 
Neglecting the full derivatives (and all possible anomalies) one can write 

Yi = ~< {qi T q) ■ {[<NsV £ q + q<J T£ [V £ q}) > 
Now, using equations of mouton Vg = it is easy to find 
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Yi = ~< (qYq) ■ ([qV T ]q - q[W T )q}) > (21) 
Finally from C-parity it is clear that 

Yi = -X, (22) 
In the same way, if we neglect anomalies, one can show 

Xi = - < (qi5l £ q) ■ ^([gV e ]7 5 g + g7s[V £ g]) >= 

= ~\< ([^]757 £ g + g757 £ [V £ )g])(g 75 <?) >= (23) 

Using Fiertz transformation both by color and scalar indexes, one can express vacuum 
averages (19) through each other. Finally a system of exact equations can be found, which 
solution, taking in account (22,23), is: 

Y 2 = -Y 2 = -X 2 = \x l] Y l = -Y l = -X 1] X 1 = X 2 = (24) 

So we see, that all vacuum averages of operators, constructed from four quarks and one 
derivative are expressed through X\. One must emphasize, that (24) is exact statement, 
don't based on factorization hypothesis. 

To estimate X ± , we use factorization hypothesis analogously to what we have done in 
sect. 1. Let us consider vacuum averages of the operator Z 1 

Z 1 =< 2{q\ n {D a G a p) n q) • (g[V p q]) > (25) 
which can be factorized as 

Z i = ~\< {q\ n {D a G a p) n [V ]) >< qq > (26) 

Now we can use equations of motion D a G^ = —gq{\ n /2)^ (for simplicity we shall limit 
us by case with one flavor, two flavor case is similar). Then, taking into account (24), we 
get: 



Z 1 = ± < {q\ n [Vpq}){q\ n ^q) X qq >= j^X 2 < qq >= ~X, < qq > -g (27) 
On the other hand, using equations of motions, (25) can be rewritten as 

Z ± = - < {q\ n q)(q^\ n q){q[V p q}) > -g (28) 

Now, using Fiertz transformation and also use C-parity and neglecting full derivatives (in 
the same way as in (21-22)), we can write 



Zl = I'll < ^)(^)^ V ^D > + < (Ql 5 l a e aPT£ q) ■ (q<y T£ q) • (flV/j?]) >} (29) 



The first term here can be factorized, so according our anzatz we have 



Z 1 = l<qq> ■X 1 -g (30) 
o 

Comparing (27) and (30) we get | X 1 — — | X 1 , So 

X 1 = (31) 

Then from (24) we can conclude, that all vacuum averages (19) are zero. One can easily 
see, that every vacuum averages of the operators with dimension 7 can be expressed 
through a set vacuum averages, discussed in this two sections. Thus, results, obtained in 
sections 1,2 allows one to give estimations for all possible vacuum averages of the operators 
with dimension 7. This estimation may be significant in a large range of problem, where 
contribution of high dimension operators became necessary. 

An example of such problem is the calculation of Bjorken [6] and Ellis -Jaffe [7] sum 
rules in the framework of QCD sum rules, offered in [8] we are going to discuss in next 
section 

Section 3 

In this section we use results obtained in the previous sections for the analysis of the 
power corrections to the first moment of the structure function of a polarized nucleon. 
On importance of the power corrections 1/Q 2 to structure functions of polarized nucleon 
was indicated in [9,10], where necessity of their contribution was considered to satisfy 
with experimental data on deep inelastic scattering on a polarized [11] nucleon, on the 
one hand, and Gerasimov - Drell - Hearn sum rule [12,13] - with other (see also discussion 
of this problem in connection with "spin-crisis" problem in review [14,15]). By the most 
natural candidate for this role seems the contribution of the operators of twist 4. The 
contribution of these operators to deep inelastic scattering on a polarized nucleon was 
calculated in [16]. 



M S(NS) 



k s ^ • Lr s \i - * s (q 2 )m - l <<uS Q 7 ]>> y 



+ 



4m|r 
3 Q 2 



/ dxx 2 (gl +n ^ n \x) + \g{ +n ^ n \x)^ + 0(l/g 4 ) (32) 



Here: 



K S(NS) = A ( Iv 

18 V' 



9l S 



G, 



1.25; g s A = 0.1 ±0.04 (see [14]) 



Gy 

< f/s^s) > are defined as < at I uS(ns) | N > = Sft < uS (ns) > 

~ 18 
U s = u G^yU + (u — > d) + — (u — > s) 
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U NS = u G^ vlv u -{u^d) 

and Sfj, = N'j^'j^N; N be a nucleoli spinor. <C JJ S ( NS ) ^> are matrix element of twist 
four, we are interest in this paper. <C U ^> has been calculated in [8] from sum rules for 
3-point correlator 

r r A 2 <^ U S( ~ NS ^ ^> 
r M (p) = i 2 / dare*" / < T[ V (x)U^ NS \y)f}(0)] >= -2p^ 5 ^— 2 ^ + - (33) 

Where X p is proton coupling and rj is proton current [17] 

V = e abc {u a C lx u b ) lblx d c 

As always in QCD sum rules, the correlator (33) is considered at large negative p 2 . How- 
ever in this case, in difference from usual 3-point correlator, though x ~ 1/p is small, but 
there are no limitations on y and, therefore, it is necessary to take into account region 
y ^> x ~ 1/p too. This lead to the fact, that except usual vacuum expectations (local 
operators) in operator expansion also appears field induced vacuum expectations - 2-point 
(bilocal) correlators of the type 



d 4 y < T{O,(y)O^0)} > 

(see papers [18 - 21], where this approach was offered and discussed). 

In some cases these bilocal operators can be reduced to local ones, using low -energy 
relations (see, for example [18,22]), in other cases one should consider corresponding 2- 
point sum rules to estimate this bilocal operators (see [20,21], for example). 

In the [8] the correlator (33) was calculated and the power correction in 1/p 2 up to 
contribution of the operators with dimension d = 8 was accounted. After borelization the 
result, obtained in [8], has the form 

« U » +R ■ M 2 = - e ^ ( ^/ M2) ^2AM 2 Jds s 2 e- s / M2 ln ^ 2 /s+ 

p V o 

+ BM\1 - exp(-s /M 2 )) + CM 2 + D^j (34) 

A, B, C, D correspond to loop contribution (A) and power correction of operators of 
dimension 4,6,8Q. (Here, and also in expressions (35, 36,40) we for simplification omit 
subscript S(NS) in all cases, where it is obvious). Note, that (34) depend on ultraviolet 
cut-off parameter fi 2 , even after borelization. This is consequence of a simple, but incorrect 
model of continuum accounting in [8] , based on ordinary dispersion relation, as was noted 
in [14] . In [14] was offered the method, how one should correctly exclude continuum, using 
double dispersion relations and was shown, that in this method dependence of unphysical 
cut-off parameter /i 2 disappear. The procedure, offered in [14] lead to following sum rules 
(instead of (34)) 

1 For the scalar case <C U s 3> contribution of s-quark in [8] was neglected. 
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« U»+R.M* = _ ™P(rn%/M>) ^ ^ ^ 



2A 2 

p y o 



2AM 2 J ds 
o 

e- s ' M \ S() + s In s /s) + M 4 |l - (1 + ^K 50 ^ + CM 2 + £> j 



(35) 



Here: 



A b = a s /n • 4/5; B b = -32/3 • tt 2 / 2 5 2 

A NS = ajTr ■ 4/9; B NS = - < ^ > 
C 5 = (Zn s /M 2 + 0.5) • 32/27 • a s /7r ■ a 2 + 8/9 • tt 2 • II 
C 7 ^ = (Zn s /M 2 + 1) • 32/27 • ajir ■ a 2 + 8/9 • tt 2 • IT 
£) s = -1/9 -ml -2 -a 2 
D NS = -l/3-m 2 -2-a 2 
a = -An 2 < jpi/j > 

R correspond to the contribution of single-pole terms, II is bilocal power correc- 
tion, which was estimated in [8] as IT = 3.10~ 3 GeV 6 ; other parameters are standard: 
U = 0.133GeV, 5 2 = 0.2 GeV 2 [23], m 2 ~ 0.8GeV 2 [17]; a s (lGeV)~ 0.37, < ifip >= 
-0.014GeV 3 , < g 2 G 2 >= 0.5GeV 2 ; A 2 = 32tt 4 A 2 = 2.1 GeV 6 , and continuum threshold 
so = 2.25GeV 2 ([17, 21], see also [14]). Hereafter we will use this correct result for <C U ^> 
(35), but one should note, that results and conclusions we will present at the end of this 
section, are similar both for (34) and (35). 

To single out « U » itself, according to [8] the operator 

d 



1 - M 2 - 



dM 2 

acting on both sides of (34) was used. Using this operator for improved result (35), we 
find 



< £/>= 



2A 2 



2A J ds s{m 2 N - s)e~ s/M \s + s In s /s) + 



- M 2 )(l - (1 + ^)e-/- 2 ) + (^) 2 e-^Y 
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+ (Cm 2 N + 32/27 • a s /n ■ a 2 ■ M 2 ) + D(l + m 2 N /M 2 ) 



(36) 



Analysis of sum rules (34), according [8], lead to < U s >~ -(0 + 0.1)GeV 2 , < 
jjns > _ o.l8 at Borel mass M 2 ~ lGeV 2 and even we take into account improved form 
(35,36), results change slightly (see [14]). However from (34-36) it is possible to see, that 
in <C U ^> just dominates the contribution of the operators of dimension 8, (both for (34) 
or (35,36)) that is the last accounted term of expansion. Thus, there is the problem on a 
reliability of results, obtained in [8], and for this purpose it is necessary to evaluate the 
following term in expansion, that is contribution of the operators of dimension D = 10. As 
usual, we shall consider that sum rules (36) are reliable, if contribution of the operators 
of dimension D = 10 will appear less than contribution of the operators of dimension 



In this section we will estimate the dimension 10 operators contribution to sum rules 
for <C U 3> (36). We will take into account only tree diagrams (fig. la-Id), because all 
other are suppressed by loop factors like 1/4-7T 2 . 

The main problem here is the estimation of vacuum averages of dimension 10 operators, 
and this can be done by help of results of sect. 1,2. In this section we'll calculate the 
contribution off all local operators of D=10 (fig.la,b) and also those bilocal operators 
(fig.lc), for which exact low-energy relations exists (see eq.(27) in [17], where this relation 
was discussed for very similar case). 

In this work we don't take in account diagrams of fig. Id, because they consist bilocal 
operator contribution, the calculation of which is a separate problem. Nevertheless the 
estimations allow to hope, that the contribution of this diagram will hardly essentially 
change an obtained result. 

It is necessary to note, that in [17] similar 3-point correlator for the current qG^^q 
was considered up to dimension 10. But in our case we take into account greater number 
of the diagrams because, using results of the previous sections and our anzatz, we can eval- 
uate practically all arising vacuum averages of dimension 10, and not just only those from 
them, which can immediately factorized to forms < g 2 G 2 >< ipip > 2 or (< ipG^a^ip >) 2 . 

Let us also do some notices about the diagram on fig. lb. On the first sight it express 
in terms of unfactorizible vacuum average like: 



where T££!f - any matrix, constructed from 7^-martrixes (and g^ v , e^ va P also). However 
up to full derivatives (37) can be written as 



D 
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K =< UxgG^VaUp] ■ U T [Vf3Ua] > -Tj 



(37) 




+ < (u\(D pG p V )\W a Up}) ■ (u r u a ) > 



+ < (u x Gp U [Vf3V a u p ]) ■ (u T u a ) > + 
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+ < (uG^[V a u p })([u T Vp}u a ) > j (38) 

Now one can easily see that all vacuum averages in right side of (37) can be factorized. 

All other diagrams (fig. la, lb) can be immediately factorized by use of results 
(18), (24), (31). Finally, result for sum rules (36), taking into account the contribution 
of operators of dimension 10 (fig.la-lc), is: 



_ e m%IMi 

< [/>= 



2A2 



so 

2A J ds s{m 2 N - s)e~ s/M \s + s In s /s) + 



+s(K-m 2 )(i-(i + ^) + (-^)' 



e so/M* + 



E ( \ m 2 N 



+ {Cm 2 N + 32/27 • a s /n ■ a 2 ■ M ) + D(l + m N / M ) + — {l + , ; ^ 2 
where contribution of dimension 10 operators are 

NS _<^> 2 ( 4 [10 2 2 11 4 \ 1 1 4.^,^2 



(39) 



E s = < ^ o >2 \m 4 +(^-< g 2 G 2 > +^rrA U < ^ > 2 ( 10) 



18 u \ 9 J 6 u / 9 



In (40) we omit all vacuum averages like Rf, or S\ or X\- < tptp > and so on, which have 
been shown in sect. 1,2 to be negligible small, and also take into account result for Rd 
from (18). 

Using standard values < g 2 G 2 >, < ijjip > one can see, that for <^ U NS ^> contribution 
of dimension 10 is only about 20-25 % of those of dimension 8 and are within the limits 
of permissible accuracy. So for <^ U NS ^> our results confirm a conclusion, made in [8], 
that contribution of operators of twist 4 in Bjorken sum rules (6) are small. 

But for case C f/ s > contribution of dimension 10 became approximately more or 
equal that contribution of dimension 8, (at M 2 = lGeV 2 ) so for <^ U s ^> it seems that 
sum rules (35,39) are inapplicable. So, most likely, for twist 4 operators contribution to 
M s in (32) (i.e. in for Ellis- Jaffe [7] sum rules) it is impossible to make any predictions 
from QCD sum rules (at least in those approach as discussed). Note also, that some other 
real reasons, that also make result for <C U s 3> doubtful are discussed in [14]. 

Author would like to thank B.L.Ioffe for very precious discussions and many useful 
advises. 
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Figure captions 



1. Diagrams, corresponding of the dim. 10 contribution. Circles denote derivatives, 
dashed lines denote gluons, solid lines denote quarks. 
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